TTP99-49 
December 1999 



Differential Equations for Two-Loop Four-Point Functions 

T. Gehrmann and E. RemiddiQ 

Institut fur Theoretische Teilchenphysik, Universitdt Karlsruhe, D-76128 Karlsruhe, Germany 



Abstract 

At variance with fully inclusive quantities, which have been computed already at the two- or 
three-loop level, most exclusive observables are still known only at one-loop, as further progress 
was hampered so far by the greater computational problems encountered in the study of multi-leg 
amplitudes beyond one loop. We show in this paper how the use of tools already employed in inclusive 
calculations can be suitably extended to the computation of loop integrals appearing in the virtual 
corrections to exclusive observables, namely two-loop four-point functions with massless propagators 
and up to one off-shell leg. We find that multi-leg integrals, in addition to integration-by-parts 
identities, obey also identities resulting from Lorentz-invariance. The combined set of these identities 
can be used to reduce the large number of integrals appearing in an actual calculation to a small 
number of master integrals. We then write down explicitly the differential equations in the external 
invariants fulfilled by these master integrals, and point out that the equations can be used as an 
efficient method of evaluating the master integrals themselves. We outline strategies for the solution 
of the differential equations, and demonstrate the application of the method on several examples. 
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1 Introduction 

Perturbative corrections to many inclusive quantities have been computed to the two- and three-loop level 
in past years. From the technical point of view, these inclusive calculations correspond to the computation 
of multi-loop two-point functions, for which many elaborate calculational tools have been developed. In 
contrast, corrections to exclusive observables, such as jet production rates, could up to now only be 
computed at the one-loop level. These calculations require the computation of multi-leg amplitudes to 
the required number of loops, which beyond the one-loop level turn out to be a calculational challenge 
obstructing further progress. Despite considerable progress made in recent times, many of the two-loop 
integrals relevant for the calculation of jet observables beyond next-to-leading order are still unknown. 
One particular class of yet unknown integrals appearing in the two-loop corrections to three jet production 
in electron-positron collisions, to two-plus-one jet production in electron-proton collisions and to vector 
boson plus jet production in proton-proton collisions are two-loop four-point functions with massless 
internal propagators and one external leg off-shell. 

Taking two-loop four-point integrals arising in the calculation of Feynman diagrams in covariant 
gauges (non-covariant gauges introduce integrals of a structure beyond the treatment of this paper) as 
an example, we elaborate on in this paper several techniques to compute multi-leg amplitudes beyond 
one loop. We demonstrate how integration-by-parts identities (already known to be a very valuable tool 
in inclusive calculations) and identities following from Lorentz-invariance (which are non-trivial only for 
integrals depending on at least two independent external momenta) can be used to reduce the large 
number of different integrals appearing in an actual calculation to a small number of master integrals. 
This reduction can be carried out mechanically (by means of a small chain of computer programs), 
without explicit reference to the actual structure of the integrals under consideration and can also be 
used for the reduction of tensor integrals beyond one loop. 

The master integrals themselves, however, can not be computed from these identities. We derive dif- 
ferential equations in the external momenta for them. Solving these differential equations, it is possible 
to compute the master integrals without explicitly carrying out any loop integration, so that this tech- 
nique appears to be a valuable alternative to conventional approaches for the computation of multi-loop 
integrals. 

The plan of the paper is as follows. In Section || we review the derivation of the integration by 
parts (IBP) identities and introduce the Lorentz invariance (LI) identities. In Section || the differential 
equations for the master amplitudes are obtained. The practical application of these tools is outlined 
in detail in Section |^ on a self-contained rederivation of the one-loop massless box integral with one 
off-shell leg. Section || contains examples of massless two-loop four-point functions with one off-shell 
leg, evaluated for arbitrary space-time dimension. We show which of these functions can be reduced to 
simpler functions and which are genuine master integrals, and compute some of the master integrals by 
solving the corresponding differential equations. Finally, Section ^ contains conclusions and an outlook on 
potential future applications of the tools developed here. The higher transcendental functions appearing 
in our results for the one-loop and two-loop integrals are summarised in an Appendix, where we also 
discuss how these functions can be expanded around the physical number of space-time dimensions. 

2 Reduction to Master Integrals 

Any scalar massless two-loop integral can be brought into the form 



where the Di are massless scalar propagators, depending on k, I and the external momenta pi, . . . ,p n while 
Si are scalar products of a loop momentum with an external momentum or of the two loop momenta. The 
topology (interconnection of propagators and external momenta) of the integral is uniquely determined by 
specifying the set (D\, . . . , D t ) of t different propagators in the graph. The integral itself is then specified 
by the powers m; of all propagators and by the selection (Si, . . . , S q ) of scalar products and their powers 




(2.1) 
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(ni, . . . , n q ) . (all the ra, are positive integers greater or equal to 1, while the are greater or equal to 0). 
Integrals of the same topology with the same dimension r = £\ m i °f the denominator and same total 
number s = J^. n, of scalar products are denoted as a class of integrals It,r,s- The integration measure 
and scalar products appearing the above expression are in Minkowskian space, with the usual causal 
prescription for all propagators. The loop integrations are carried out for arbitrary space-time dimension 
d, which acts as a regulator for divergencies appearing due to the ultraviolet or infrared behaviour of the 
integrand (dimensional regularisation, ^ ||). 

Any four-point function depends on three linearly independent external momenta, p±, P2 and ^3. 
At the two-loop level, one can combine the two loop momenta k and I and these external momenta to 
form 9 different scalar products involving k or I. As the propagators present in the graph are (linearly 
independent) combinations of scalar products, only 9 — t different scalar products can appear explicitly 
in an integral with t different propagators. Since a two-loop four-point function can have at most seven 
different propagators, as can be found by considering the insertion of a propagator into a one-loop 
four-point function, one has in general t < 7, while the minimum number of massless propagators in 
a two-loop graph is t = 3, corresponding to a two-point function. The number of different two-loop 
four-point integrals for given t (number of different propagators), r (sum of powers of all propagators) 
and s (sum of powers of all scalar products) can be computed from simple combinatorics: 

m,r,.) = ( r - 1 M s - ,+s )- w> 



The number N(I t , r ,s) of the integrals grows quickly as r, s increase, but the integrals are related among 
themselves by various identities. One class of identities follows from the fact that the integral over the 
total derivative with respect to any loop momentum vanishes in dimensional regularisation 

d d k d , , . 

(^^-)=0- (2-3) 

where J is any combination of propagators, scalar products and loop momentum vectors. J can be a 
vector or tensor of any rank. The resulting identities || |j| are called integration-by-parts (IBP) identities 
and can for two-loop integrals be cast into the form 

d d k d d l d , , 

«"/(*:, I, Pi) = 0, 



(2w) d (2ir) d dk» 

d d k d d l d 
{2ir) d (2^81^ 



u"/(M,Pi) = 0, (2.4) 



where the integrand f(k,l,pi) is a scalar function, containing propagators and scalar products and 
can be any external or loop momentum vector. As a consequence, one obtains for a graph with m loops 
and n independent external momenta a total number of A^ibp = m(n + m). For a two-loop four-point 
function, this results in ten IBP identities for each integrand. 

In addition to the IBP identities, one can also exploit the fact that all integrals under consideration 
are Lorentz scalars (or, perhaps more precisely, "d-rotational" scalars) , which are invariant under a 
Lorentz (or (i-rotational) transformation of the external momenta. In order to derive the resulting Lorentz 
invariance (LI) identities, we consider an infinitesimal Lorentz transformation 

p" -> p" + 5p" =p fl + 5^p v with Set = -5e^ , (2.5) 

which should not change the scalar Feynman integral 

Z(pi + Spx, .. . ,p n + Sp n ) = J(pi, . . . ,p n ) . (2.6) 

Expanding 

d d 

J(pi + 5pi,...,p n + Sp n ) = Z(pi, . . . ,p n ) + 5p^-^I(pi, . . . ,p n ) + ... + 5p^-^I{pi, ■ ■ ■ ,p n ) , (2.7) 
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one arrives at 

d „ d 



P ^ + --- +P -Wn) IiPl '-- Pn) ^ - (2 - 8) 

Since 5e^ has six independent components, the above equation contains up to six LI identities. These are 
however not always linearly independent. To determine the maximum number of linearly independent 
identities, one uses the antisymmetry of Se^ to obtain 

opifi dpi,, dp nil opnuj 

This equation can be contracted with all possible antisymmetric combinations of Pi^Pj V to yield LI 
identities for I. For a three-point vertex, two of the external momenta are linearly independent (n = 
2), and only one antisymmetric combination of them can be constructed, resulting in one LI identity 
{Njji = 1). A four-point function depends on three external momenta (n = 3), allowing us to construct 
three linearly independent antisymmetric combinations, which yield three LI identities (iVu = 3). The 
full potential of the LI identities can only be exploited for integrals involving five or more external legs, 
which allow to construct six linearly independent antisymmetric combinations of external momenta, thus 
projecting out all six LI identities (Ali = 6). 

Since / is a scalar, it can not depend on the momenta pi itself, but only on scalar products Sij = 2pi -pj 
of the external momenta. Replacing 

£2^+^)^, (2.10) 



one finds that (2.9) becomes a trivial identity, independent of I. However, the derivatives in (2.9) can 
be interchanged with the loop integrations in /, s uch that they do not act anymore on the integral I, 
but on the integrand of /. After this interchange, ( [2.9] ) becomes a non-trivial relation between different 
integrals. 

In the case of two-loop four-point functions, one has a total of 13 equations (10 IBP + 3 LI) for each 
integrand corresponding to an integral of class It, r ,s, relating integrals of the same topology with up to 
s + 1 scalar products and r + 1 denominators, plus integrals of simpler topologies (i.e. with a smaller 
number of different denominators). The 13 identities obtained starting from an integral It,r,s do contain 
integrals of the following types: 

— It, r ,s- the integral itself. 

— It-i, r ,s'- simpler topology. 

— It.r+i,s, It.r+i,s+i ■ same topology, more complicated than It, r ,s- 

— 7{ )r _i )S , It t i — same topology simpler than It,r,s- 

Quite in general, single identities of the above kind can be used to obtain the reduction of /t, r +i,s+i 
or I t , r +i,s integrals in terms of It, r ,s and simpler integrals - rather than to get information on the It,r,s 
themselves. 

Integrations-by-parts identities are widely applied in multi-loop calculations of inclusive quantities 
(see e.g. Q for a review), which are related to two-point functions. In these calculations, only a relatively 
small number of different topologies has to be considered, but the integrals appearing in the calculation 
can bear large powers of propagators and scalar products, arising for example from expansions in masses 
or momenta. In these calculations, it is desirable to have reduction formulae for arbitrary powers of 
propagators and scalar products. These can be obtained from IBP identities derived for an integral with 
arbitrary powers (to be treated symbolically) of scalar products and propagators; the derivation of these 
symbolic reduction formulae requires a lot of ingenuity, based on the direct inspection of the explicit form 
of the IBP identities for each considered topology, and could not be carried out mechanically. 

For loop integrals with a large number of external legs, IBP identities are needed for a large number of 
different topologies, but in general for relatively small powers of propagators and scalar products. In this 
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case, it would therefore be desirable to have a mechanical procedure for solving, for any given topology, 
IBP and LI identities for integrals with fixed powers of the propagators and scalar products. 

If one considers the set of all the identities obtained starting from the integrand of all the N(It, r , s ) in- 
tegrals of class It, r ,s, one obtains (iViBP -\-Nu)N(It, r ,s) identities which contain N(It r +i a +\)+N(It r +i a ) 



integrals of more complicated structure. From (2.2) it can be read off that with increasing r and s the 
number of identities grows faster than the number of new unknown integrals^. As a consequence, if 
for a given i-topology one considers the set of all the possible equations obtained by considering all the 
integrands up to certain values r* , s* of r, s, for large enough r*, s* the resulting system of equations is 
overconstrained and can be used for expressing the more complicated integrals, with greater values of 
r, s in terms of simpler ones, with smaller values of r, s. (Let us observe that, the system being overcon- 
strained, the equations cannot be all independent; it is not a priori known how many equations are in 
fact linearly independent and, correspondingly, how many integrals of the topology under consideration 
will remain after reduction). 

The required values r* and s* for r and s can be found by counting the number of accumulated 
equations (equations for all integrals with r < r* and s < s*) and comparing them with the number of 
accumulated unknown integrals, with (r < r* + 1, s < s* + 1), but excluding (r = t,s = s* + 1). As 
an example, we list in Table ^ the number of equations and unknowns for two-loop four-point functions 
with seven denominators t = 7. It can be seen that a complete reduction requires at least one of the 
combinations (r*,s*): (7,2); (8,1); (9,0). 



t = 7 



different I, 



t,r,s 



s 

r 
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2 


3 
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14 


21 


28 


35 
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28 


56 


84 


112 


140 


10 


84 


168 


252 


336 


420 



, , , equations 
accumulated 

unknowns 



\^ s 
r 





1 


2 


3 


7 


13 
22 


39 
45 


78 
76 


130 
115 


8 


104 
106 


312 
213 


624 
354 


1040 
535 


9 


468 
358 


1404 
717 


2808 
1196 


4680 
1795 



Table 1: Comparison of the number of (IBP and LI) identities to the number of new unknowns (different 
integrals with t = 7) appearing in these equations for a two-loop box integral with t = 7 internal 
propagators. The identities for It, r ,s contain at most It.r+i,s+i- It can be seen that for growing r and 
s, the number of equations, upper number in each box, exceeds the number of unknowns, given by the 
lower number. 



The above table illustrates that typically hundreds of equations have to be solved in order to obtain 
a reduction towards simpler integrals. The task is performed automatically (and independently of the 
topology!) by a computer program invoking repeatedly the computer algebra packages FORM || and 
MAPLE §. For any given four-point two-loop topology, this procedure can result either in a reduction 
towards a small number (typically one or two) of integrals of the topology under consideration and 
integrals of simpler topology (less different denominators), or even in a complete reduction of all integrals 
of the topology under consideration towards integrals with simpler topology. Left-over integrals of the 
topology under consideration are called irreducible master integrals or just master integrals. If a topology 
turns out to contain irreducible master integrals, one is in principle free to choose which integrals are 
taken as master integrals, as long as the chosen integrals are not related by the IBP and LI identities. In 
our reduction, we choose It,t,o f° r topologies with one master integral and It,t,o together with the required 
number of integrals of type It,t+i,o for topologies with more than one master integral. 

As a final point, it is worth noting that the procedure described above can be used to reduce a 

1 The importance of this fact was first pointed out by S. Laporta and exploited in [B|. 
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tensor integral towards scalar integrals. All integ rals appearing in the projection of an arbitrary tensorial 
integral onto a tensor basis will take the form (|2.l| ), i.e. they can be classified as It,r,s and reduced to 
master integrals. 

In the following section, we will demonstrate that the established reduction of all the integrals to a 
few master integrals can be used also to write differential equations in the external invariants for the 
master integrals themselves, and then how the differential equations can be used to compute these master 
integrals. 



3 Differential Equations for Master Integrals 



The IBP and LI identities discussed in the previous section allow us to express integrals of the form (2.1) 
as a linear combination of a few master integrals, i.e. integrals which are not further reducible, but have 
to be computed by some different method. 

For the case of massless two-loop four-point functions, several techniques have been proposed in the 
literature, such as for example the application of a Mellin-Barnes transformation to all propagators [||, || 
or the negative dimension approach ||lo| . Both techniques rely on an explicit integration over the loop 
momenta, with differences mainly in the representation used for the propagators. So far, these techniques 
were only applied to a limited number of master integrals: Smirnov |^| has recently used the Mellin-Barnes 
method to compute the planar double box integral for the case of all external legs on shell (massless case) ; 
the same method has been applied by Tausk || for the computation of the non-planar on-shell double 
box integral; the negative dimension approach has been applied by Anastasiou, Glover and Oleari |jfO| to 
compute the class of two-loop box integrals which correspond to a one-loop bubble insertion in one of the 
propagators of the one-loop box. A general method for the computation of the master integrals appearing 
in two-loop four-point functions has however not yet been found. So far, it has also not even been clear 
(apart from the planar double box topology, where Smirnov and Veretin have recently demonstrated 
the reducibility of any integral of this topology to two master integrals JlT|| ) how many master integrals 
exist for a given topology. Solving the identities discussed in the previous section, we are now able to 
identify the irreducible master integrals. A list of reducible two-loop four-point topologies will be given 
in Section |a. 

A method for the computation of master integrals avoiding the explicit integration over the loop 
momenta is to derive differential equations in internal propagator masses or in external momenta for the 
master integral, and to solve these with appropriate boundary conditions. This method has first been 



suggested by Kotikov |12 to relate loop integrals with internal masses to massless loop integrals. It 
has been elaborated in detail and generalised to differential equations in external momenta in Mj; first 
applications were presented in [ fl4| . In the case of four-point functions with one external off-shell leg and 
no internal masses, one has three independent invariants, resulting in three differential equations. 

The derivatives in the invariants Sij = (pi + Pj) 2 can be expressed by derivatives in the external 



momenta: 



d 1 / ^ d M d ^5 



Sl2 a^ = 2\ + ^dpJ +P2 M~ P *dp» 



ds 23 2 V dp 1 ^ dp% dp% 

The combinations of derivatives and momenta appearing on the right hand side of ( |3.f| ) are obviously 
linearly independent from the combinations appearing in the LI identities ( |2.9| ) (which vanish identically 
when acting on a function depending on the scalars S12, s%3, S23). The three derivatives of an integral 
It,r,s(si2, S13, S23, d) are not linearly independent, but related due to the properties of I under rcscaling 
of all external momenta: 

It,r,.(«12, *13, S23, d) = h,r,s (A 2 Sl2 , A 2 S 13 , A 2 S 23 , d) , (3.2) 
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where a{d, r, s) is the mass dimension of the integral. For a m-loop integral in d space-time dimensions 
with r powers of denominators and s scalar products, one finds a(d,r,s) = md + 2s — 2r. The above 
equation yields the rescaling relation 



a 
2 



Sl2 



d 
ds\2 



Sl3 



d 



S23 



d 



ds 



23 



It,r,s(si2, S13, S23, d) = 



(3.3) 



Let us now first consider the case of a topology with only one master integral, which is chosen to be 
h,t,o, defined as 

f d d k d d l 

It,t,o( s ij> s jk,Ski,d) = J —-^——f tjtfi (k,l,Pi) ' ^'^ 

where ft,t,o{k, l,Pi) is a suitable integrand of the form appearing in ( |2.1| ). It is evident that acting with the 
right hand sides of (3.1) on the right hand side of (3.4) will, after interchange of derivative and integration, 
yield a a combination of integrals of the same type as appearing in the IBP and LI identities for It,t.o, 
including integrals of type It,t+i,i and It,t+\,o- Consequently, the scalar derivatives (on left hand side of 
(3.1)) of It,t,o can be expressed by a linear combination of integrals up to It,t+i,i and It,t+i,o- These can 
all be reduced to it,t,o and to integrals of simpler topology by applying the IBP and LI identities. This 
reduction results in differential equations for It,to of the form: 







Sij q — I t> t,o(sij,Sjk, s ki , d) = A(sij , Sjk, s^, d)I ttt fi(sij , Sjk, su, d) 

OSij 



+ F(s lj ,Sjk,S k i,d, I t -l, r , s (Sij , Sjk, Ski, d)) 



(3.5) 



where sy, Sjk, Ski are the three invariants, A{sij,Sjk,Ski,d) turns out to be a rational function of the 
invariants and of d, F(It-\, r ,s) is a linear combination (with coefficients depending on Sy, Sjk, Sfcj and 
d) of integrals of type It-i, r , s , containing only topologies simpler than It,t,o, but potentially with high 
powers of some denominators and scalar products. These I t -i^ : s, which refer to a simpler topology and 
therefore can be considered as known in a bottom-up approac h, p lay the role of an inhomogeneous term 
in the equation; one can then look for the proper solution of (|3.5| ) in a straightforward way. 

The master integral It,t,o{sij, Sjk, Ski, d) can indeed be obtained by matching the general solution of 



(3.5) to an appropriate boundary condition. Quite in general, finding a boundary condition is a simpler 
problem than evaluating the whole integral, since it depends on a smaller number of kinematical variables. 
In some cases, the boundary condition can even be determined from the differential equation itself: for 



s^ = 0, ( p.5| ) yields, if ,4(0, s jk ,s k i,d) ^ 0, 



It,t, o(0, sjk, Ski, d) = - [A(0, Sjk, s^, d)] 1 ^(0, Sjk,Ski,d,I t -i, r , s (0, Sjk,su,d)) 



(3.6) 



For A vanishing at s%j — one can consider It t t,o m t ne limit where one of the external momenta vanishes, 
corresponding to the vanishing of both invariants involving this momentum, e.g. s^j = and s k i = for 
p% = 0. In this case, It,t,o reduces to a three-point vertex function with one off-shell external leg. All these 
functions, which might be determined by iterating the procedure just described for the considered 4-point 
function, have actually already been computed at the two-loop level in |l5[ using IBP identities to reduce 
all possible topologies to a few master integrals, which, in this case, can be computed straightforwardly 
using Feynman parameters (cf. Section ^|). Starting from the boundary condition in = Ski = 0, one 
can determine It,t,o{®, s jk, Ski,d) by solving the differential equation in Ski - this provides the desired 
boundary condition in Sy = 0. 



For topologies with more than one master integral, (3.5) will be replaced by a system of coupled, 
linear, first order differential equations for all master integrals of the topology under consideration. The 
determination of the master integrals from these equations follows the same lines as discussed above, with 
the only difference that the general solution for the system of coupled differential equations is harder to 
obtain than for a single equation. In case the coupled equations can not be decoupled by an appropriate 
choice of variables, several mathematical techniques can be employed here fil| : the system of n coupled 
first order differential equations can for example be rewritten into one n-th order differential equation, 
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which is then solved with standard methods. In some cases, the system can also be transformed into a 
form which is known to be solved by generalised hypergeometric series jl6], [l7|, [l8|, |l9| . 

It is clear from the above discussion, that the determination of a master integral of a certain topology 
with t different denominators requires the knowledge of all the integrals appearing in the inhomogeneous 
term. These integrals are subtopologies of the topology of the integral under consideration, and contain at 
most t — 1 different denominators. The determination of master integrals has therefore to proceed bottom- 
up from simpler topologies with a small number of different denominators towards more complicated 
topologies with an increasing number of different denominators. For the case under special consideration, 
massless two-loop four-point functions with up to one external leg off-shell, this implies that one has 
to progress from the simplest master integrals with t = 3 (off-shell two-point function) to construct all 
master integrals up to t = 7. 



4 A Pedagogical Example: the One-Loop Four-Point Function 

To illustrate how the method explained above works in practice, we present in this section a detailed and 
self-contained derivation of the one-loop four-point function 



q 




P2 


pi 




P3 



d d k 



1 



(2ir) d k 2 (k - p 2 ) 2 (k - p 2 -P3) 2 {k -px~P2 ~ P3) 2 



(4.1) 



The topology of this integral is given by the set of its four propagators, it is the only t = 4 topology at 
one loop. In the notation introduced above, (4.1) corresponds to an integral 14,4,0, the simplest integral of 
this topology with all propagators appearing in first power and with no scalar products. For the one-loop 
four-point function, one has four independent scalar products involving k^ (k 2 and A^»pf ) and four linear 
independent denominators. This implies that any scalar product involving k^ can be rewritten as linear 
combination of propagators and invariants Sij = 2pi M p^. Consequently, no integrals with scalar products 
in the numerator can appear for this topology. 

The reduction of all four integrals with one squared propagator (14,5,0) can be carried out by consid- 
ering the IBP identities for 14,4,0: 



d d k d 



(2ir) d dk» k 2 {k - p 2 ) 2 {k -p 2 - p 3 ) 2 {k -px-p2- P3) 2 



= 0, 



(4.2) 



where can be the loop momentum k^ or any of the external momenta pf, thus yielding four identities. 
It turns out that LI identities and IBP identities for integrals with higher powers of the propagators do 
not contain additional information which would allow a reduction of 14,4,0- The integral (4.1) is therefore 
a master integral. It is the only master integral for this topology. 

To proceed towards the differential equations in the invariants s\ 2 , s\ 3 and s 23 for this master integral, 
let us consider the derivatives in the external momenta: 



Pi 



d 
dpt 



I- 1 

Pa 



dp% 



q 




P2 


pi 




P3 



q 




P2 


pi 




P3 



q 




P2 


pi 




P3 



A d k 



1 
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/ 2pi(fc - Pi - Pi zBtk 
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& d k 



(4.3) 



(27r) d k 2 {k -p 2 ) 2 (k -p 2 -P3) 2 (k-pi - pa - P3) 2 
{2p 2 l (k-p 1 -p 2 -p 3 ) f _ l 2p 2 l {k -p 2 -pa)n . 2p 2 I (fc-p 2 ) 



V (k-Pi 
d d k 



P2 -P3) 2 



(k - p 2 - p 3 ) 2 
1 



(k-P2? 



,(4-4) 



(27r) d k 2 (k -p 2 ) 2 {k - p 2 -p 3 ) 2 (k, -p x -p 2 -p 3 ) 2 
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(fc ~pi - p 2 - p 3 ) 2 (k-p2-P3) 2 J 

The right hand sides of the above equations contain terms with 

(i) four different propagators with one squared propagator, no scalar product 

(ii) four different propagators, no squared propagator, no scalar product 

(iii) three different propagators, one squared propagator. 

The terms of type (i) can now be reduced to type (ii) and (iii) by using the integration-by-parts identities 
(4.2). One obtains: 
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where (•) denotes a squared propagator and (x) stands for a pinched (cancelled) propagator. 
A check on these equations is provided by the rescaling relation (3.3) 
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which is related to the above equations by 
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Inserting the derivatives obtained above, one finds that the rescaling relation is indeed fulfilled. 

The three-propagator terms appearing in the above equations can be further reduced by using the 
IBP identities for the corresponding vertex; the vertex amplitudes, in turn, can all be expressed in terms 
of bubble integrals, and one finally finds 
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A similar identity is obtained by exchanging p\ «-> p2- 

The differential equations for the one- loop box integral then follow from (3.1). The set of differential 
equations reads: 
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where S123 = S12 + S13 + S23. The one-loop bubble diagrams in the inhomogenous term yield: 
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The boundary conditions in Sy = can be readily read off from the above: 
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. (4-14) 
. (4.15) 



The result for the one-loop box integral can in principle be obtained by integrating any of the differential 
equations ( |4.9| )-(4.11). In practice, it turns out to be more appropriate to introduce a new set of variables, 
namely S13, S23 and S123 = S12 + S13 + S23, corresponding to the arguments appearing in the two-point 
functions in the inhomogeneous terms. This transformation yields a differential equation in S123, which 
will be used for integration. Note that this transformation also modifies the differential equations in S13 
and S23- The differential equation in S123 reads: 
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(4.16) 



The boundary condition in S123 = can not be trivially determined from this equation, reflecting the 
fact that the massless one-loop box integral with all external legs on shell is not reducible to simpler 
subtopologies by IBP identities. The boundary condition in S123 = — S13 — S23 can however be determined 
from ( gig ). 



Equation (4.16) is a linear, inhomogeneous first order differential equation of the form 
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9 



which can be solved by introducing an integrating factor (see for instance p0| or any standard book on 
differential equations) 

M(x)=eS f( - x)dx , 

such that y(x) = l/M(x) solves the homogenous differential equation (g(x) = 0). This yields the general 
solution of the inhomogenous equation as 
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where the integration constant C can be adjusted to match the boundary conditions. 
For (Lit), we have at once the integrating factor 
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This factor is not unambiguous, since 
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would also be a valid integrating factor. We select ( 4.17 ) by requiring a real integrating factor in the 
region: — s 123 > — s 13 — s 23 > 0. T he final result for the box integral does not depend on the selection of 
the integrating factor; using ( 4.18[ ), one must however be more careful in applying analytic continuation 
formulae and in multiplying non- integer powers of the invariants. 

With the integrating factor (4.17), the one- loop box integral reads: 
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From a computational point of view, the evaluation of the box amplitude of (|4.lj) has been reduced to 
the one dimensional integration corresponding to the solution of (4.16). 

The lower boundary of the integral is independent of S123 and can be adjusted arbitrarily. The first 
two terms in (4.19) can be easily integrated by shifting the integration variable to s' 123 — S13 — s 2 3, which is 
then integrated between and S123 — S13 — S23. To integrate the last term, one introduces a new variable 
s i23( s 'i23 — s i3 — s 23)j which is integrated between and Si 2 3(si23 — S13 — S23). The resulting integrals 
yield can be identified as integral representation of the hypergeometric function 2-Fi- With this choice 
of variables and boundaries, no constant term is required to match the boundary conditions. The result 
for the one-loop box integral then reads: 
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The invariants appearing in this expression can be safely continued from the region — S123 > — S13 — S23 > 
to the physical region S123 > S13 + S23 > 0. The arguments of the hypergeometric functions are ratios 
of invariants, they are not changed by the analytic continuation. The non-integer powers of invariants 
appearing as coefficients acquire imaginary parts, their signs are uniquely determined by the convention 
— p 2 = —p 2 — iO, thus fixing the imaginary part of the whole expression. The above equation reproduces 
the well-known result from the literature, e.g. plf . It should be kept in mind that in applying dimensional 
regularisation, no distinction between infrared and ultraviolet poles, which both show up as l/(d — 4) 
in the above equation and in the coefficient ^2,1,01 is made. Needless to say, those singularities are an 
intrinsic feature of the dimensional regularisation, and by non means an artifact of the differential equation 
approach. In particular, the most singular part of the above one loop integral, which is ultraviolet finite, 
as well as of the two loop integrals discussed in the following section, arises from soft configurations, and 
can in principle be re-derived by applying a strong ordering procedure [p2[ to the integrands. 

When the box integral is expressed in the above form, where no expansion around d = 4 has yet been 
performed, analytic continuations, e.g. to the on-shell case S123 = or to collinear and soft limits Sy = 0, 
can be made with ease. 



5 Results on Two-Loop Four-Point Functions 

In the following, we shall outline how the techniques derived in Sections ^ and |^ can be applied to 
the computation of two-loop integrals appearing in amplitudes for the decay of one massive into three 
massless particles: two-loop four-point functions with one off-shell leg. The main purpose of this section 
is to illustrate applications of the tools developed above to non-trivial problems; the list of integrals given 
here is far from complete. We provide a comprehensive list of master integrals only up to t — 5, for t = 6 
and t = 7 only reducible integrals are quoted. 

A prefactor common to all massless scalar integrals is 
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T(d - 3) 
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which is also appearing in all counterterms in the MS-scheme. 

In the following, the notation of external momenta is as follows: pi denotes an on-shell momentum, 
Pij(k) denotes an off-shell momentum, being the sum of two (three) on-shell momenta pi, pj(, pk) with 
ij(k) = (Pij{k)) 2 - P is an arbitrary momentum. 



5.1 t = 3 

For two-loop integrals with t = 3, only one topology exists: the two-loop vacuum bubble. The corre- 
sponding integral fulfils a homogeneous differential equation, which can not be used to infer any boundary 
condition. The integral can however be computed using Feynman parameters: 
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Several different two-loop topologies exist for t = 4. Two types of two-point functions are encountered. 
The first can be reduced to 
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(5.3) 



using IBP identities. The second is the product of two one- loop bubble integrals and yields 
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which generalises trivially to a three-point function 
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Only one of the master integrals at t = 4 fulfils a homogeneous differential equation: a two-loop vertex 
integral with one off-shell leg. This integral can also be computed using Feynman parameters: 
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The other vertex integral topology with one off-shell leg can be reduced using IBP identities: 
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Among the three vertex integrals with two off-shell legs, only one can be reduced using IBP and LI 
identities: 
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The two remaining ones are master integrals. Written out in terms of propagators, they read: 
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Both fulfil inhomogeneous differential equations. For a vertex P123 —> P12 +P3, the appropriate variables 
for the differential equations are S123 and s\ 2 . To illustrate the structure of the differential equations, we 
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(5.11) 



The boundary conditions for S123 = or S12 = are obtained directly from the vertex integrals with one 
off-shell leg quoted above. Using these, one finds 
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The second master integral can be obtained from this by analytic continuation of the hypergeometric 
function: 
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Vertex integrals with three off-shell legs can not appear as subtopologies in two-loop four-point func- 
tions with one off-shell leg. 



5.3 t = 5 

The two-loop two-point function with t = 5 is a well known example || for the application of IBP 
identities: 
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The four different t = 5 three-point functions with one off-shell leg can also be reduced by using IBP 
and LI identities: 
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By applying IBP and LI identities, it is likewise possible to reduce all but one t = 5 three-point 
function with two off-shell legs: 
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The remaining three-point function is a master integral, which can be found by solving the corre- 
sponding differential equations: 
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For t = 5, one finds four different topologies for two-loop four-point functions. These are all master 
integrals obeying inhomogeneous differential equations in the external invariants. Solving these equations 
and matching the boundary conditions, all master integrals can be determined. 
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The integrals ( 5.25 ) and ( p.2q ) are one- loop bubble insertions into the one-lo op b ox and hav e already 
been computed for arbitrary powers of the propagators in jn^. The integrals ( 5.27 ) and ( [3.28] ) were, to 
our knowledge, not known up to now. In the reduction of integrals of the topology (5.28), one finds two 
master integrals, whose differential equations decouple in the variable A = S13 — S23. The second master 
integral for this topology can be found by rearranging one of the differential equations: 
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Finally, products of one- loop vertex with one- loop bubble integrals also yield topologies with t = 5. 
These can all be reduced to ( ^5 ) and are not quoted explicitly. 

The complete list of integrals at t = 5 which were derived in this section can now be used to compute 
all integrals at t = 6 and t = 7 which can be reduced using IBP and LI identities. The results of this 
reduction are summarised in the following. 



5.4 t = 6 

Two-loop integrals with t = 6 arising in calculations in covariant gauges must be three- or four-point 
functions. Since we are concerned with subgraphs that can appear in the reduction of four-point functions 
with one off-shell leg, we need to consider three-point functions with up to two off-shell legs. For general 
three-point functions at t = 6, one finds three distinct topologies: two planar and one crossed arrangement 
of the loop momenta. The crossed graphs correspond to master integrals, while the planar graphs are 
reducible, as first pointed out in Jl5|| , where the three-point integral with one off-shell leg was computed. 
We reproduce these results: 
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(d-3)(3d-8)(3d-10) 1 P23 

(d-4) 3 S 23 (S12 + S13)S123 

4(d- 3) (3d- 8) (3d- 10) 1 ,} 




P123 




-P23 
-Pi 



(d - 4) 3 
3(d-3)(3d- 10) 



(S12 + si 3 )s 2 2 3 
P123, 





+ 



+ 



2(d-4)(d-5) (si2 + s 13 ) 2 

3(d- 3) (3d- 8) (3d- 10) 1 

2(d-4) 2 (d-5) S2 3 (si2 + s 13 ) 2 

3(d- 3) (3d- 8) (3d- 10) 1 



-P23 
-Pi 



P123, 




P23 




Pl23, 




(5.30) 



(5.31) 



P^rr 



(d-4) 2 (Sl2 + S1 3 )S123 

-*—Pi 



-Pi 

-P23 



(5.32) 
-P23 
-Pi 



2(d-4) 2 (d-5) si2 3 (si2 + si 3 ) 2 
For t — 6, both one-loop bubble insertions on propagators of the one-loop box can be reduced: 



(5.33) 



P123- 
Pi- 



-P2 

-P3 



3(si2 + S13) 
S13S23 



P123- 



Pi- 



+ 



3(3d-10) 1 
2(d-4) S13S23 



P23 




-P2 

-P3 
-P2 



P123, 



-P3 




-P2 
-Pl3 



17 



3(3d-8)(3d- 10) 1 



2(d-4)2 s 2 3fi23 




Pl23- 



-P2 



-P3 



3 
Sl3 



Pl23- 



-P2 



+ 



+ 



3(3d- 8) (3d - 10) 



-P3 

1 /'in 



3(3d - 10) 1 
2(d-4) S13S23 




(5.34) 
-Pi 
-P23 




2(d-4) 2 S13S23S123 

Among the remaining four planar diagrams at t = 6, three are reducible to simpler subtopologies: 



(5.35) 



P123- 



P123- 



Pi- 



P123- 



Pi- 



-P2 



-P3 



-P2 



-P3 



-Pi 



-P3 



Pl23- 



= -3- 



512 
S13S23 



Pi-*- 



-Pi 
-P3 



6(d-3) 1 
d-4 si 3 



P123- 



Pi- 



-Pi 

-P3 



3(d- 3) (3d- 10) 1 P123, 



(d-4)* 



S13S23 




-Pi 



-P23 



3(d- 3) (3d - 8) (3d - 10) 



(d-4) 3 



Pl3 



S 13 S23 




P23 



S13S| 3 




Pl23- 



= -3- 



£12 
S13S23 



Pi- 



-Pi 

-P3 



+ 



6(d - 3) £12 + S13 
d-4 S13S23 



P123- 



Pi- 



3(d-3)(3d-10) 1 P23 



(d - 4) 2 



S13S23 




+ 



3(d-3)(3d-10) 1 P!23, 



(d-4)2 



S13S23 




-P2 

-P3 
-P2 

-P3 
-P2 

-Pl3 



3(d-3)(3d-8)(3d- 10) 



(d - 4) 3 



P23 



S13S23 



P123- 




Pl23, 



Sl3S23 s 123 




d-4 



Pi- 



-Pi 
-P3 



(5.36) 



,(5.37) 
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+ 



6(d - 3) £12 + S13 
d-4 S13S23 



P123- 



Pi- 



P23 




+ 



3(d- 3)(3d- 10) 1 
(d-4) 2 S13S23 

3(d-3)(3d-10) 1 Piaa 



-P2 
-P2 



(d - 4) 2 



S13S23 




-»— P3 
-»— P2 



-P13 



3(d- 3) (3d - 8) (3d - 10) 



P13 



(d-4) 3 



S 13 S 23 




(5.38) 



One of the two remaining non-planar diagrams is also reducible, the other non-planar topology contains 
two master integrals. The reducible integral reads: 



P123- 
Pi- 




-Pi 
-Ps 



+ 



+ 



+ 



+ 



+ 



3(d- 


4) 


S12 


2d- 


9 


S13S23 


3(d- 


4) 


S23 


2d- 


9 


S12S13 


3(d- 


4) 


«13 



P123- 

Pl- 
P123- 



P123-*- 



-Pl 

-P3 
-P2 



2d - 9 S12S23 



Pi- 



-*~Pi 

-*~P3 



-P2 



3(d - 3) (3d - 8) (3d - 10) s 13 + s 23 Pl2 



(d-4) 2 (2d-9) 



Sl2 s 13S23 




3(d-3)(3d-8)(3d- 10) s 12 + s 23 Pl3 

(d-4) 2 (2d-9) s 12 s 2 3S23 

3(d - 3) (3d - 8) (3d - 10) s 12 + s 13 P23 

(d-4) 2 (2d-9) si2Si3S 2 3 

3(d-3)(3d-8)(3d-10) 1 P123, 





(d-4) 2 (2d-9) S12S13S23 




(5.39) 



5.5 t = 7 



At i = 7 different one finds six different topologies. Three of them are triangle insertions to the one-loop 
box. These three integrals are all reducible, two of them contain only master integrals up to t = 5 in 
their reduction: 



P123- 

Pi- 



-P2 

-P3 



6(d-3)(3d-14) 1 
(d-4)(d-6) 4~3 



P123- 
P2- 



-Pi 
-P3 
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+ 



3(3d-14) s\ 2 
(d - 6) s 2 13 4 3 



Pl23- 



Pi- 



-Pi 
-P3 



6(d-3)(3d-14 ) ( Sl2 + s 13 ) 2 
(d-4)(d-6) 



Pl23- 



„2 ,2 
A 13 A 23 



Pi-*- 



-Pi 



-*~P3 




+ 



3(d - 3) (3d - 10) (3d - 14) (2d - 10)si 2 + (3d - 14)s 23 Vii3 ( 
2(d-4) 2 (d-5)(d-6) s 13S 2 3 ( Sl2 + S23 ) 

3(rf-3)(3rf-10)(3d-14) (2d_10)s 12 + (3d - 14)s 13 Viz 
2(d-4) 2 (d-5)(d-6) 



2 „2 



13 a 23 




-P2 

-P13 
-P2 

-P3 



3(d - 3)(3d - 10)(3d - 14) (2d - 10)si 2 .si 2 3 + (3d - 14)s 13 s 23 Pi*} 
2(d~4) 2 (d-5)(d-6) 5f 3S 2 3 ( Sl2 + S23 ) 

3(d - 3) (3d - 8) (3d - 10) (3d - 14) (2d - 10)si 2 + (d - 6)s 23 Pl3 




-*—P2 

> 

-*—Piz 



2(d-4)3(d-5)(d-6) 



"l3 S 23( S 12 + s 23) 




+ 
+ 

P123- 



3(d - 3) (3d - 8) (3d - 10) (3d - 14) (d - 5)s 12 - (d - 4)s 13 P23 




(d-4)3(d-5)(d-6) sf 3 ^ 3 
3(d - 3)(3d - 8)(3d - 10)(3d - 14) (2d - 10)si 2 s 123 + (d - 6)si 3 s 23 Pl23 , 



Pi- 



2(d-4) 3 (d-5)(d-6) 

-P2 



S 13 S 23 S 123(S12 + S23) 




-P3 



2(2d-9) 2s 12 + s 13 + s 23 



P123- 



(d-4)(d-6) Si 3 S 23 

P123- 



Pl- 



-P2 



-P3 



3(d-4) ( S i 3 + s 23 ) 2 



d- 6 



„2 2 
S 13 S 23 



Pl- 



-P2 



-P3 



6(d-3)(3d-14) (S12 + S23) 2 



13 a 23 



(d-4)(d-6) 

6(d-3)(3d- 14) (S12 + S13) 2 
(d-4)(d-6) S ? 3S 2 3 



P123- 

P2- 
P123- 

Pl- 



-Pl 

-P3 
-P2 

-P3 



+ 



3(rf-3)(3rf- 10)(3ri- 14) (2d - 10)s 12 + (3d - 14)s 23 Pw f 
2(d-4)2(d-5)(d-6) S 2 3S 2 3 



-Pi 

-P3 



, (5-40) 
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3(d - 3) (3d - 10) (3d - 14) (2d - 10)si 2 + (3d - 14)si 3 P^j 
2(d-4) 2 (d-5)(d-6) 

3(d- 3)(3d- 10) (3d - 14)si 2 + (4d - 18)si 3 - (d - 4) s 23 ^ 
' (d-4)2(d-6) s^sg 

3(d - 3) (3d - 10) (3d - 14)s 12 - (d - 4)si 3 + (4d - 18)s 23 
' (d-4)2(d-6) 



-P2 
-P3 




-P13 




-»— i?2 

P23 



-Pi 



3(d - 3)(3d - 8)(3d - 10) (d - 5)(3d - 14)(s 12 + s 13 ) + (d - 4) 2 s 23 Pia 
(d-4)3(d-5)(d-6) ^ 

3(d - 3)(3d - 8)(3d - 10) (d - 5)(3d - 14)(s 12 + s 23 ) + (d - 4) 2 s 13 P23 
(d-4) 3 (d-5)(d-6) 




„2 3 
*13 5 23 




(5.41) 



The remaining three topologies are the double box and two different momentum arrangements of the 
crossed box. These topologies contain each two master integrals. 



6 Conclusions and Outlook 

Progress in the computation of exclusive observables, such as for example jet production rates, beyond 
the next-to-leading order has up to now been hampered mainly by difficulties in the calculation of virtual 
two-loop integrals with more than two external legs. In contrast to this, many inclusive observables 
(which correspond from the calculational point of view to two-point functions) are known to next-to- 
next-to-leading order and even beyond. These higher order calculations relied on a variety of elaborate 
technical tools for the computation of the virtual integrals. In this paper, we outline how techniques 
known from multi-loop calculations of two-point integrals can be modified and extended towards the 
computation of integrals with a larger number of external legs. 

We demonstrate that the large number of different two-loop integrals appearing in an actual calculation 
can be reduced to a small number of scalar master integrals by using the well-known integration-by-parts 
identities |j| together with identities following from Lorentz-invariance which are unique to multi-leg 
integrals. As a by-product of this reduction, one is also able to reduce two-loop integrals with tensorial 
structure to scalar integrals. In contrast to two-point integrals, where only a few topologically different 
graphs can appear with potentially large powers of propagators and scalar products, one finds that the 
reduction of three- and four-point integrals gives rise to a large number of topologically different graphs, 
which appear however only with small powers of propagators and scalar products. The reduction of two- 
point functions usually proceeds via solving manually the integration-by-parts identities for arbitrary 
powers of propagators and denominators in a given graph topology; this procedure seems to be not 
practicable for multi-leg integrals. To accomplish the reduction of these, we developed an algebraic 
program which automatically derives and solves the integration-by-parts and Lorentz-invariance identities 
for a given graph up to some pre-selected fixed number of powers in denominators and scalar products 
independent of the topology. 

To compute the scalar master integrals, we derive differential equations in the external momenta Jj"3"[ 
for them; the boundary conditions of these differential equations correspond to simpler integrals, where for 
example one of the external momenta vanishes. These differential equations can be solved (for arbitrary 
space-time dimensions) by employing standard mathematical methods. We observe that the differential 
equations for the master integrals we considered up to now are solved by generalised hypergeometric 
functions. We illustrate the application of this method in detail on the example of the one-loop four- 
point function with one off-shell leg. 

Using the differential equation method, we provide a complete list of all master integrals with up 
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to t — 5 different denominators that can appear in the reduction of two-loop four-point functions with 
one off-shell leg. We also list all reducible integrals with t = 6 and t = 7 different propagators. The 
computation of the master integrals with t = 6 and t = 7 is still an outstanding task. 

The differential equations for these outstanding master integrals are of similar structure as the differ- 
ential equations for master integrals with a smaller number of different propagators. The main problem 
towards a complete computation of these integrals is at present the integration of the inhomogeneous 
term, containing itself already hypergeometric functions arising from the subtopologies. 

It is worthwhile to point out similarities and differences between the differential equation method 
employed in this paper and other methods employed for similar calculations in the literature. Both the 
negative dimension approach of |ic(| and the Mellin-Barnes transformation method employed in Ji], ^) rely 
on choosing a particular assignment of momentum vectors to the internal loop propagators. After this 
assignment, a representation of the propagators in terms of a multiple sum (negative dimension approach) 
or an integral transformation (Mellin-Barnes method) is employed, such that the integral over the loop 
momentum can be carried out explicitly. The final result for the integral is then retrieved by resummation 
of a multiple sum or by an inverse integral transformation. Both methods, when employed for arbitrary 
space-time dimension, give rise to generalised hypergeometric functions, which can be represented as 
multiple sums as well as inverse Mellin-Barnes integrals Jl7| . In the differential equation method, one 
assigns momentum vectors to the loop propagators only for the sake of deriving the differential equations 
and the IBP and LI identities. After applying these identities to simplify the differential equations, 
one obtains a relation between the derivative (with respect to an external momentum) of a master 
integral, the master integral itself and other master integrals with simpler topology, independent of 
the parametrisation chosen for the internal propagators. These differential equations can then be solved 
analytically by integration; the resulting integrals correspond to the integral representations of generalised 
hypergeometric functions |]l7| . Using the differential equation method, one can therefore circumvent the 
explicit loop momentum integration needed in the other methods and one arrives at a representation of 
the hypergeometric functions, which is presumably more transparent than a multiple sum or an inverse 
integral transformation. In the integral representation, it is in particular straightforward to identify 
linear combinations of different hypergeometric functions, which are difficult to disentangle in the other 
representations. At present, it should however not be claimed that any of the methods is superior, since 
none of them could yet be employed to compute all outstanding two-loop four-point master integrals. 

As a final point, we note that the methods derived in this paper contain a high level of redundancy, 
which allows for a number of non-trivial checks on the results obtained with them. The automatic 
reduction to master integrals using integration-by-parts and Lorentz-invariance identities corresponds to 
the solution of a linear system of equations containing more identities than unknowns. The existence of 
a solution to this system provides therefore already a check on the self-consistency of the identities. In 
computing the master integrals from differential equations, one integrates one of the three differential 
equations in the external invariants, such that the result can be checked by inserting it in the remaining 
two differential equations. 

In short, this paper demonstrates how techniques developed for multi-loop calculation of two-point 
functions can be extended towards integrals with a larger number of external legs. As a first example 
of the application of these tools in practice, we computed some up to now unknown two-loop four-point 
functions, relevant for jet calculus beyond the next-to-leading order. The most important potential 
application of these tools is the yet outstanding derivation of two-loop virtual corrections to exclusive 
quantities, such as jet observables. 
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A Special Functions 

This appendix summarises the series and integral representations of the hypergeometric functions appear- 
ing in the master integrals. The properties of these functions, in particular their region of analyticity, 
their analytic continuation as well as reduction formulae, can be found in the literature [Tt], [j^, 
Hypergeometric functions are sums with coefficients formed from Pochhammer symbols 

, . T(a + n) . . 

(a) n = (A.1) 

The hypergeometric functions of a single variable are given by: 

2 F 1 (a, b- c-z) = ^ — — , (A.2) 

n— 

TP I h h ^ (a)n(bl)n(b2)n z n 

3 F 2 {a,b 1 ,b 2 ;c 1 ,c 2 ;z) = > , . . r r . (A.3) 

^ ci)„ c 2 )„ n! 



Two types of hypergeometric functions of two variables also appear in our results 

(a)m+n(bl) m (b 2 )n Z^ Z. 
—_ ,, (c) m +n 

(ffll)m+n(l2)m+n(&)m Zj" 2» 



-Pi ( a i bi,b 2 ;c; z±,Z2) = > 7 t r-r , (A.4) 

— [C) r ~ 1 - m ' n ' 



5i (a 1 ,a2,b;c,d;zi,z 2 ) = F~\ 7J\ FT ' ( 5 ) 

m.n— 



These functions have the following integral representations: 

m r 

T(b)T(c-b)J 



Fi(a,b;c;z) = ^...^ / dt t 6_1 (l - 1) 6-6 -^! - te) _ ° 



3 F 2 (a,6i,6 2 ;ci,c 2 ;2) 



Re(6) > 0, Re(c - b) > (A.6) 

r(ci)r( C2 ) 



r(&i)r(d - &i)r(fc)r(c2 - fc) 

dti di 2 - ti) !-^- 1 ^ - ta) *- 6 "- 1 ^ - ht 2 z)~ a 

Jo 

Re(6i) > 0, Re(ci - 61) > 0, Re(6 2 ) > 0, Re(c 2 - b 2 ) > (A.7) 



Si (ai,o 2 ,6;c, d; z l7 z 2 ) 



T(a)T(c - a) J 

Re(a) > 0, Re(c - a) > (A.8) 

r(ai)r(c - ai)r(6)r(d - 6) 

f dtj / d^f "^(l - ti) - 01 - 1 ^ - t,)^ 6 - 1 (1 - iiMi - ti2 2 )- a2 
Jo Jo 

Re(ai)>0, Re(c-ai)>0, Re(6) > 0, Rc(d - 6) > (A.9) 



B Expansion of Hypergeometric Functions 

To separate divergent and finite parts of the loop integrals derived in this paper, one has to expand them 
around the physical number of space-time dimensions in the parameter e = (4 — d)/2. We demonstrate 
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in this appendix, that this expansion can, at least for the hypergeometric functions in one variable, 
be carried out in a mechanical way, giving rise to harmonic polylogarithms (HPL), a generalisation of 
Nielsen's polylogarithms |2^] introduced in |p4| . 

Expanding the integral representation of 2-P1 in e yields simple powers of (t, 1 — t, 1 — tz) times the 
product of some number of (hit, ln(l — t), ln(l — tz)). The powers of (t, 1 — £, 1 — tz) can be integrated 
by parts until one obtains non-trivial integrals 

d*(i, T ^ ir L ; )ln-*ln-.(l-t)In«.(l- to ). 

All these integrals are combinations of harmonic polylogarithms H(a;z), where a is a vector of indices 
with w = rii + n 2 + n 3 + 1 components, w is called the weight of the harmonic polylogartihm. The proof 
by induction in w is trivial, once the HPL formalism p4[ is recalled: 

1. Definition of the three HPLs at w = 1: 

H(l;z) = -ln(l-z), 

H(0;z) = Ins, 

H(-l;z) = ln(l + z) (B.l) 



and the three fractions 



f(0;z) 

f(-i;z) 



1 



1 - 
1 

? 

z 

1 



1 + z 



such that 



2. For w > 1: 



which results in 



— H(a;z) = f(a;z) with a = +1,0,-1 



ff(O,...,O;20 = ^ln^z, 

H(a,b;z) = / dxf(a;x)H(b;x), 



-^H(a,b;z) = f(a;z)H(b;z) 



B.2) 
B.3) 

B.4) 
B.5) 

B.6) 



This last relation is a convenient tool for verifying identities among different HPLs. Such identities 
can be verified by first checking a special point (typically z = 0) and subsequently checking the 
derivatives. If agreement in the derivatives is not obvious, this procedure can be repeated until one 
arrives at relations involving only HPLs with w = 1. 

3. The HPLs fulfil an algebra (see Section 3 of p4|), such that a product of two HPLs (with weights w\ 
and W2) of the same argument z is a combination of HPLs of argument z with weight w — wi + W2- 

Using these properties of the HPL, one can show that the integrals appearing in the e-expansion of 
the hypergeometric function can be reexpressed as 

(B.7) 
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Following the argumentation of Section 7 of j24[ , one can show that the integral on the right hand side 
of the above equation yields a linear combination of HPLs of weight w = w a + + 1. The proof goes 
via induction in Wb- 

For w a — w\, — one has H(a; t)H(b; zt) = 1. The i-integral in (B/7) yie lds then a combination 
of HPL of weight w = 1 (HI). Likewise, for Wb = the right hand side of (B/7) will yield a linear 
combination of HPLs of weight w = w a + 1 and of argument z, as proven in Section 7 of p4| . 

Considering 

f d<[-,— !— , 1 , | H(a;t)H(B,b;zt) , (B.8) 



we observe that 



(9- 



d dt 
dzt-l/z 

^H(B,b;zt) 
oz 



t 1 -t t- l/z 



1 

^2 t" 



— dt — — h boundary terms 

-l/z dt 



tf(B;zt)H(b; zt) . 



Making these replacements in ( B.8 ) and applying partial fractioning to all denominators, we are left with 



1 1 



1 



z 1 — z 1 + z 



dt - 



1 1 



1 



t 1 -t t- l/z 



H(a; t)H(b; zt) 



(B.9) 



which is a combination of HPL s with argument z and weight w = w a + Wb + 1 multiplied with (l/z, 1/(1 — 
z), 1/(1 + z)). Integrating ( B.8 ) over z will thus yield a combination of HPLs with argument z and weight 
to + 1, which completes the proof by induction. 

The e-expansion of 3F2, corresponding to a double integral in t\ and £2, is obtained by carrying out 
the procedure described here twice, again resulting in a combination of HPLs. A systematic e-expansion 
of F\ and Si, which are functions of two variables z\ and z^, will in general go beyond the harmonic 
poly logarithms in one variable. 
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